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Abstract 

The emission of hard real photons from thermalized expanding hadronic matter is dominated 
by the initial high-temperature expansion phase. Therefore, a measurement of photon emission in 
ultrarelativistic heavy-ion collisions provides valuable insights into the early conditions realized in 
such a collision. In particular, the initial temperature of the expanding fireball or equivalently the 
equilibration time of the strongly interacting matter are of great interest. An accurate determination 
of these quantities could help to answer the question whether or not partonic matter (the quark gluon 
plasma) is created in such collisions. In this work, we investigate the emission of real photons using 
a model which is based on the thermodynamics of QCD matter and which has been shown to 
reproduce a large variety of other observables. With the fireball evolution fixed beforehand, we are 
able to extract limits for the equilibration time by a comparison with photon emission data measured 
by WA98. 



1 Introduction 

In the hot and dense system created in an ultrarelativistic heavy- ion collision (URHIC), the relevant 
momentum scales for typical processes taking place inside the strongly interacting matter drop as a 
function of proper time r: Initially, the relevant scale is set by the incident beam momentum, leading 
to hard scattering processes which presumably can be described by perturbative Quantum Chromody- 
namics (pQCD). Secondary inelastic scattering processes subsequently lower the momentum scales due 
to particle production. At later times, equilibration sets in and typical momenta p are determined by 
the temperature T of a given volume clement as (p) = 3T. As the matter expands, energy stored in 
random motion of particles (temperature) is transferred to collective motion (flow), leading to a de- 
screase of T with r. Therefore, by selecting an observable associated with a given momentum scale, one 
simultaneously selects a time period in the evolution of the system. 

For this reason, a measurement of hard real photons is an ideal tool to study the early moments of the 
fireball expansion: High momentum photons are not only sensitive to early proper times, but, being 
subject to electromagnetic interactions only, their mean free path in the fireball matter is also much 
larger than the spatial dimension of the system, therefore they are capable of leaving the emission region 
without significant rescattering. Therefore, hard photons complement a measurement of low mass, low 
momentum dileptons which are dominantly emitted near the kinetic freeze-out point. 

Ideally one would like to use hard photon emission as a thermometer to determine the initial temper- 
ature reached in an URHIC and use this information to verify the creation of a quark-gluon plasma 
(QGP). However, in reality one has to disentangle thermal contributions to the photon spectrum from 
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contributions coming from initial hard scattering processes. The interpretation of the photon spectrum 
alone can therefore not be unambiguous. 

In this work, we compare a model calculation of photon emission from a fireball created in an 158 
AGeV Pb-Pb collision at SPS with data obtained by the WA98 collaboration [J. In a recent paper 
[2], we have developed a fireball model which is based on information from hadronic observables and 
lattice QCD thermodynamics, as manifest in a quasiparticle picture of the QGP. This model has been 
shown to successfully describe low mass dilepton emission 2 and, within the framework of statistical 
hadronization, the measured abundancies of hadron species [3]. In the present work, we demonstrate 
that the same model is also capable of describing the observed photon emission. The fixed setup of the 
model also enables us to establish constraints on the equilibration time To, which entered the model on 
an ad-hoc basis so far. 

This paper is organized as follows: First, we introduce the photon emission rate used in the calculation 
and discuss its interpretation in the framework of the quasiparticle picture of the QGP which has been 
used in the fireball evolution model. In the next section, we summarize the main properties of the 
evolution model and discuss constraints for the initial expansion phase. Afterwards we present the 
resulting photon spectrum and demonstrate that, in agreement with our expectation, hard photons 
originate dominantly from the early evolution phase. We investigate the possibility of using the photon 
data to set limits on the equilibration time of the fireball matter and conclude by comparing to the 
results obtained by other groups. 



2 The thermal photon emission rate 



2.1 The QGP contribution 



As we expect the dominant contribution to the spectrum of hard photons to come from a region of high 
temperatures, we focus on the QGP rate. To leading order, this rate can be calculated by evaluating 
the four diagrams shown in Fig^ Here the last two diagrams are promoted to leading order because of 
near-collinear singularities. 




Figure 1: Leading order processes for photon production in the QGP (from left to right): QCD Compton 
scattering, qq annihilation, Bremsstrahlung and annihilation with scattering (aws). 



The complete calculation of the rate to order a s has been a very involved task which has been finished 
only recently 0] |SJ H3 |S1 El E3 • F° r the present calculation, we use the parametrization of the rate 
given in [TJ]|. There, the rate for photons of momentum k is written as 

"' A 1 :A(k) (ln(T/m q (T)) + ± hx{2E/T) + C tot (E/T)) , (1) 
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with E = k and m 2 (T) = 4ira s T 2 /3 the leading order large momentum limit of the thermal quark mass. 
The leading log coefficient A(k) reads 



A(k) = 2aN c Y,<I 2 s^f 1 fD(E). (2) 

s 

Here, the sum runs over active quark flavours and q s denotes the fractional quark charges in units of 
elementary charge. The Fermi-Dirac distribution Jd{E) dominates the momentum dependence of the 
rate: To a good approximation it decreases exponentially with E. The dependence on the specific photon 
production process is contained in the term Ctot(E /T): 

C tot (E/T) = C 2 ^ 2 (E/T) + C brems {E/T) + C aws (E/T). (3) 

All these functions C(E/T) involve non-trivial multidimensional integrals which can only be solved 
numerically. In |1U) . parametrizations for the results are given as 

C 2 ^2(E/T) ~ 0.041(.E/T) -1 - 0.3615 + 1.01 exp[-1.35E/T] (4) 

and 

(E/T)+C aws (E/T) ~ 

/[ 1~~~ / 0.548 ln(12. 28 + 1/{E/T)) 0.133S/T \ (5) 

V + 6 f [ {EjrW 2 + VI + (E/T)/W27 J ' 



2.2 A quasiparticle interpretation 

In 0, we have used a picture of massive, non- interacting quark and gluon quasiparticles to describe the 
QGP. Close to the phase transition these quasiparticles are subject to confinement, parametrized by a 
universal function C (T) which reduces the number of thermodynamically active degrees of freedom as 
n{T) = n {T)C(T) with n Q (T) = J ^ f B[D) {E p /T). Here, f B(D) {E p /T) denotes the Bose (Fermi) 
distribution. In |11|. it has been shown that this ansatz is capable of describing the lattice results for 
the QCD thermodynamics extremely well. 

In the case of photon emission, we cannot strictly retain this interpretation: The process qq — > 7 is 
kinematically impossible and all other emission processes would be suppressed by a em for non-interacting 
quasiparticles. On the other hand, the fact that we want to study hard photons implies that at least 
one of the particles in the initial state is also hard. But such particles penetrate the screening cloud of 
thermal fluctuations which is ultimately responsible for the notion of weakly interacting quasiparticles. 
Therefore it appears reasonable to allow for interactions of plasma particles with momenta well above 
the scale set by the temperature. 

The remaining properties of the quasiparticle approach are encoded in the quasiparticle mass m(T) 
and the 'confinement factor' C(T). In the limit of large temperatures, the quasiparticle mass in |11| is 
chosen such as to coincide with the result of hard thermal loop (HTL) resummed calculations for the 
thermal self-energy. Therefore, inserting massive quarks as degrees of freedom into the above results 
would amount to double counting, since those already incorporate HTL resummation, at least as long 
as we consider only temperatures above 1.5 Tc- 

There is no equivalent of the confinement factor C(T) in the calculations described in 01 El El El El E|- 

We can estimate the effect of introducing C(T) as follows: 

A typical diagram, say qq annihilation, which contributes Ro to the total emission rate has the structure 
Rq ~ Jd(E /T) 2 \M\ 2 (1 + /b(E/T)), with the thermal quark distributions fn in front of the squared 
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matrix element M. corresponding to the process in vacuum and a Bose enhancement factor for the gluon 
emitted into the final state. The modification of the rate R with respect to the rate Rq in the presence 
of C{T) will read R ~ C*(T) 2 f D (E /T) 2 \M\ 2 {1 + C{T)f B {E/T)) which is always larger than C(T) 3 R . 
In the case of quarks in the final state, C(T) leads even to a reduced Pauli blocking (the final state 
modification becomes (1 — fr>(E/T)C(T)) which is larger in the presence of confinement). 

On the other hand, as mentioned above, at least one of the incoming particles has to have a large 
momentum. Such high-momentum particles are not subject to confinement along with the bulk of the 
matter. One would expect them to hadronize well outside the thermalized region, leading to C(l) = I 
at one of the incoming legs in the diagram. 

Therefore, C(T) 3 can be taken as a conservative estimate of the effect of the confinement factor on 
the emission rate. If it can be shown that hard photon emission is dominated by a region where the 
temperature is so large that C(T) as 1, the above expression for the emission rate can be used to 
approximately describe photon emission from a system of quasiparticles also. This is also the region 
where we expect the mass of quasiparticles to be given by the HTL result. We will verify this property 
a posteriori. 

Clearly, the prescription outlined here has to be regarded as an approximation till a more detailed version 
of a quasiparticle description of the QGP incorporating confinement is available. 



2.3 The hadronic contribution 

As the temperatures in the hadronic evolution phase of the fireball are lower than in the QGP phase, 
we expect the hadronic contributions to the emission of hard photons to be small. Therefore, we will 
not discuss this contribution in great detail. 

Vector mesons play an important role for the emission of photons from a hot hadronic gas. The first 
calculation of such processes has been performed in |12| in the framework of an effective Lagrangian. It 
has been found that the dominant processes are pion annihilation, 7r + 7r - — > fry, 'Compton scattering', 
tt^p — > tt^j and p decay, p — > 7r + 7r~7. 

Several more refined approaches have been made since then (for an overview, see In the following, 

we will use a parametrization of the rate from a hot hadronic gas taken from |14| which is given as 

E J N [fin' 4 GcV- 2 ] ~ 4.8T 2 - 15 exp[-l/(1.35£T) a77 ]expb£/T]. (6) 



2.4 The integrated rate 

In order to compare to the experimentally measured photon spectrum jj, we have to integrate Eq. 
over the space-time evolution of the fireball, 

= -f- / drR 2 (r) / dz / *,nr- (7) 

In this expression, R(t) stands for the radial expansion of the fireball, Ay denotes the rapidity interval 
covered by the detector, y(z) is the rapidity of a volume element at position z and the limits of the 
k z integration come from the fact that a photon emitted at the (boosted) edge of the fireball has to 
have a longitudinal momentum in a certain range in order to be detected in the rapidity window of the 
experiment. In this expression, we have assumed spatial homogeneity and a cylindrical fireball. 



dN 
d 2 k t dy 
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3 The fireball evolution model 



In this section, we briefly outline the general framework of the fireball evolution model. The model is 
described in greater detail in [2] and [T5|. 

The underlying assumption of the model is that the strongly interacting matter produced in the initial 
collision reaches thermal equilibrium at a timescale tq w 1 fm/c. For simplicity, we assume spatial 
homogeneity of all thermodynamic parameters throughout a 3-volume at a given proper time. The 
evolution dynamics is then modelled by calculating the thermodynamic response to a volume expansion 
that is consistent with measured hadronic momentum spectra at freeze-out. 

The volume itself is taken to be cylindrically symmetric around the beam (z-)axis. In order to account for 
collective flow, we boost individual volume elements inside the fireball volume with velocities depending 
on their position. As flow velocities in longitudinal direction turn out to be close to the speed of light, 
we have to include the effects of time dilatation. On the other hand, we can neglect the additional 
time dilatation caused by transverse motion, since typically v± <C v z . The thermodynamically relevant 
volume is then given by the collection of volume elements corresponding to the same proper time r. In 
order to characterize the volume expansion within the given framework, we need first of all the expansion 
velocity in longitudinal direction as it appears in the center of mass frame. For the position of the front 
of the cylinder, we make the ansatz 

z(t)=v t + c z fdt' (8) 

where c z is a free parameter. The time t starts running at to such that z§ = vq to is the initial longitudinal 
extension with vq the initial longitudinal expansion velocity. The longitudinal position z(t) and t itself 
define a proper time curve r = \Jt 2 — z 2 (t). Solving for t = i(r) one can construct z(r) = z(t). Then 
the position of the fireball front z(t) in the center of mass frame can be translated into the longitudinal 
extension L(t) of the cylinder on the curve of constant proper time r. One obtains 

At the same proper time we can define the transverse flow velocity and construct the transverse extension 
of the cylinder as a circle of radius 

R(t) =R Q + C1 _ fdr' fdr" ^Q- (10) 

J T J To £ l T I 

where Rq corresponds to the initial overlap radius, and cj_ is a free parameter. With the values of L and 
R obtained at a given proper time we can parametrize the 3-dimensional volume as 

V(t) = vri? 2 (r) L(t). (11) 

In principle the model is fully constrained by fixing the freeze-out time tf (at the fireball front), the 
proper time ry and the two constants c± and c z , such that the measured hadronic observables are 
reproduced. In practice this is achieved only at SPS energy, where the freeze-out analysis [TQ allows a 
complete characterization on the basis of hadronic dN/dy and mt-spectra and HBT radii. For the 5% 
most central Pb-Pb collisions one requires 

R{t } ) = R Si v ± (T f ) = v{, v z {r f )=vl and T( Tf )=T f (12) 

while making a trial ansatz for the ratio p/e. Note that the proportionality of the acceleration to this 
ratio (which is reminescent of the behaviour of the speed of sound) allows a soft point in the EoS to 
influence and delay the volume expansion. 
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Assuming entropy conservation during the expansion phase, we fix the entropy per baryon from the 
number of produced particles per unit rapidity. Calculating the number of participant baryons as 

N p {b) = J d 2 s n p (b,s) = J d 2 s T A {s) [1 - exp(-a%T B (\b - s|)] + (T A «-» T B ) , (13) 

we find the total entropy Sq. The entropy density at a given proper time is then determined by s — 
Sq/V(t). Using the EoS given by the quasiparticle approach JT], thereby providing the link with lattice 
QCD, we find T(s(t)) and also p(r) and e(r), which in the next iteration replace the trial ansatz in the 
volume parametrization. Finally we arrive at a thermodynamically self-consistent model for the fireball 
which is by construction able to describe the hadronic momentum spectra at freeze-out. 

In order to find the fireball evolution relevant for the photon emission measured for the 10% most 
central collisions, we follow the procedure outlined in |2j. In brief, we consider an effective system which 
starts out with a reduced number of participants and hence reduced total entropy content. Neglecting 
azimuthal asymmetries, we keep parametrizing the expanding system as a cylinder with reduced initial 
radius. Assuming that the freeze-out temperature is approximately unchanged for more peripheral 
collisions, we determine a reduced proper evolution time and modify the geometrical freeze-out radius 
and the transverse flow velocity accordingly. However, going from the 5% to the 10% most central 
collisions, we find differences in the early evolution phases on the level of a few percent only. 

Thus, the fireball evolution is completely constrained by hadronic observables. In 2 , it has been shown 
that this scenario is able to describe the measured spectrum of low mass dileptons, and in [3] it has been 
demonstrated that under the assumption of statistical hadronization at the phase transition temperature 
T<7, the measured multiplicities of hadron species can be reproduced. None of these quantities is, however, 
sensitive to the detailed choice of the equilibration time tq. Therefore, we have only considered the 
'canonical' choice To — 1 fm/c so far. The calculation of photon emission within the present framework 
provides the opportunity to test this assumption and to limit the choice of tq. 



4 Results 

The result of the evaluation of Eq. 0) with the fireball evolution model described in the previous section 
is shown in Fig. [5] 
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Figure 2: Thermal photon spectrum for 10% most central Pb-Pb collisions at SPS, 158 AGeV Pb-Pb 
collisions, shown are calculated rate (total, contribution from QGP and hadronic gas) and experimental 
data p. 
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The overall agreement with the data is remarkably good. Above 2 GeV, the calculation underestimates 
the data somewhat, leaving room for a contribution of prompt photons from initial hard processes of 
about the same magnitude as the thermal yield. Note that the spectrum is almost completely saturated 
by the QGP contribution — for k t > 3 GeV, the hadronic contribution is almost two orders of magnitude 
down. This can in essence be traced back to the strong temperature dependence of the emission rate 
normalization and justifies the approximate treatment of the hadronic contribution a posteriori. 

In order to study the importance of the initial, high temperature phase in more detail, we present the 
time evolution of the spectrum in Fig. [3] 
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Figure 3: The total photon emission spectrum and the integrated rate at proper times r = 0.2,0.5, 1.0 
and 2 fm/c of the fireball evolution. 

One observes that the large k t region is almost exclusively dominated by the first fm/c of evolution proper 
time, whereas the yield in the low k t region is not yet saturated after 2 fm/c. However, the temperature 
associated with these evolution times is always larger than 250 MeV, leading to C(T) > 0.9 m] and 
C(T) 3 » 0.7. This justifies neglecting the effect of quasiparticlc properties on the rate a posteriori: 
There is a 30% uncertainty introduced into the calculation, but this is comparable with other intrinsic 
uncertainties, such as the detailed choice of the numerical value of as(T) in this temperature regime. 

The high kt tail of the spectrum is potentially capable of providing information about the initial tempera- 
ture reached immediately after equilibration. This capability is seriously limited in practice, however, by 
the need to assess an unknown contribution of prompt photons, which may be large in this region. Bear- 
ing this uncertainty in mind, we can nevertheless pursue this idea further in Fig. 21 where we investigate 
the sensitivity of the result to the equilibration time To of the fireball. 

We find that the low k t region of the spectrum is hardly affected by different choices for the equilibration 
time, while for larger fc t one is increasingly sensitive to short evolution timescales. An equilibration time 
of 0.5 fm/c corresponding to an initial temperature of 370 MeV leads to a good description of the data 
without the inclusion of any prompt photon contribution. On the other hand, a rather slow equilibration 
corresponding to to = 2 fm/c and an initial temperature of 260 MeV requires a sizeable contribution 
from prompt photons. 

Without any reference to prompt photons, we are therefore able to fix r = 0.5 fm/c as the lower bound 
for the equilibration time: Shorter timescales would lead to thermal photon emission overshooting the 
data. 

If we want to find an upper limit for the equilibration time, we have to address the issue of prompt photon 
emission. For this purpose, we use the results of two different works ^| E| to illustrate the possible 
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Figure 4: The thermal photon emission spectrum for different choices of the equilibration time tq as 
compared to experimental data pQ. 



range of predictions dependent on the average value of pt, an 'intrinsic' transverse momentum scale 
which is introduced as a phenomenological parameter to account for non-perturbative effects. (Fig. EJ. 
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Figure 5: Prompt photon production in Pb-Pb collisions as a function of the photon transverse mo- 
mentum kt for different values of average parton intrinsic transverse momentum (pf) |17j as compared 
to experimental data £Q (figure adapted from |17jV 

One observes that prompt photon production is able to explain the data above 3 GeV for (pf) ~ 1.8 
GeV 2 , but no choice of (p 2 ) leads to a description of the data below 3 GeV. Thus, there is a momentum 
region between 2 and 2.5 GeV where the data are only weakly affected by a prompt photon contribution. 

On the other hand, the present calculation of thermal photon emission indicates sensitivity to the 
equilibration time in this region (see Fig.^J. One finds that an equilibration time of 2 fm/c is unable to 
describe the data even in the presence of a sizeable prompt photon contribution, hence we may regard 
this as the upper limit for the equilibration time given the results of |17j . 

If we use alternatively the results of |18j to estimate the contribution of prompt photons, the situation 
is somewhat less clear, though the main conclusions remain valid. In the momentum region near 2 GeV, 
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the result obtained for the prompt photon contribution in |18| is about 50% of the thermal emission 
spectrum obtained with tq = 1 fm/c in the present calculation, which in turn is already on the lower end 
of the error bars of the data. Therefore, the prompt photon contribution alone is unable to explain the 
spectrum in this momentum region, and a thermal contribution with To ~ 2 fm/c is needed to describe 
the data near fc t = 2 GeV. (In [T^], it is stated that the prompt photon contribution is of the same size 
as the thermal emission spectrum corresponding to a fireball with initial temperature of 300 MeV, which 
is not true in the present calculation. However, this statement is derived using Bjorken hydrodynamics, 
whereas the present model incorporates longitudinal acceleration of the fireball matter which leads to 
less cooling initially and hence to a prolonged QGP phase.) 

However, there are several uncertainties in the calculation. We have already discussed the complications 
introduced by the quasiparticle picture and the detailed choice of as. An additional potential problem 
in the calculation is that the rate has been calculated to order as, but as is not a small quantity in 
the relevant temperatur region, so there might be sizeable higher order corrections. In view of these 
uncertainties and the less urgent need for an addition al contribution to prompt photons if the results 
of ^Hl are taken, it is possibly better to estimate 0.5 fm/c < To < 3 fm/c for the equilibration time. 

In principle, one might try to construct a scenario with a large equilibration time, using the most 
optimistic estimate for the prompt photon contribution. Below temperatures ~ 200 MeV, the approxi- 
mations made to calculate the photon emission from the QGP phase break down, however qualitatively 
we expect the confinement factor C(T) to strongly reduce the emission rate. Likely candidates for photon 
yield between 2 and 2.5 GeV are therefore only the hadronic evolution phase and the pre-equilibrium 
phase. 

In order for the hadronic phase to contribute significantly, either the emission rate or the four-volume of 
emitting matter needs to be increased significantly. The four-volume of fireball matter in the hadronic 
phase, however, is tightly constrained by the measured amount of dilepton radiation, which has been 
discussed in the present framework in 0. An increase of the emissivity of a hadronic gas, on the other 
hand, would very likely be accompanied by a change in the number of active degrees of freedom, which 
presumably are driven by in-medium mass reductions of resonances as argued in |19| . A strong increase 
of active degrees of freedom, however, is not in agreement with a statistical hadronization analysis done 
in 3 . Even an overall mass reduction of 10% (excluding the pseudoscalars) is not compatible with the 
observed hadron ratios. Therefore, within the present framework, the hadronic phase is not likely to 
give a large contribution to the photon spectrum. 

It remains the question of the yield from the pre-equilibrium phase. No rigorous scenario leading to 
thermal equilibrium for SPS conditions has been developed so far, however, several aspects of the pre- 
equilibrium dynamics have already been investigated. 

In |20|. the kinetic equilibration of different quark flavours was investigated under the assumption that 
gluons come to an early equilibrium and constitute a heat bath in which quark motion takes place. 
An early non-equilibrium distribution of quarks would of course directly influence the photon emission 
spectrum. Somewhat related is the question of chemical equilibration of quarks. Here, the hot-glue 
scenario [22 nas been suggested where an initial undersaturation of the quark densities with respect to 
the thermal equilibrium densities is assumed, i.e. almost all of the entropy of the system is carried by 
the gluons, leading to a drastically increased initial temperature. 

The findings of (20j suggest, albeit for RHIC conditions, that the typical timescale for the kinetic equili- 
bration of light quark flavours is of order ~ 1 fm/c. This timescale is roughly in line with the assumptions 
made in the present work. However, it leaves the question if the photon emission signal is affected if one 
starts with a suitable out-of-equilibrium initial quark distribution. 

Regarding the hot-glue scenario, note that the drastically increased temperature of the partonic matter 
would mostly affect the high momentum tail of observed photons and therefore leave the momentum 



region below 2.5 GeV less affected. This is especially true since also the overall normalization of the 
emission rate is reduced with respect to equilibrium conditions due to the undersaturation of the quark 
densities. 

A different approach to pre-equilibrium dynamics has been taken in |22j . Here, transverse momentum 
dependence of dilepton emission has been calculated in a kinetic framework and calculations for different 
initital parton distributions have been tested. While the approach is very interesting, it is hard to directly 
estimate its possible influence on photon emission within the present framework. 

In |23|. an investigation of non-equilibrium photon emission has been carried out using a parton-cascade 
model (PCM). Here, the essential findings were that only a very dilute partonic medium is created in 
the collision. Photon emission from this medium was shown to explain the data above 3 GeV when 
integrated up to the hadronization point, but below 3 GeV, the photon spectrum from the partonic 
phase falls below the data. It is difficult to relate these findings directly to the present approach, since 
no equilibrium phase in either partonic or hadronic phase is described in the PCM. However, we may take 
this as an indication that pre-equilibrium dynamics is most likely to strongly affect the high momentum 
region of the photon spectrum only where we find considerable uncertainties with regard to the question 
of intrinsic pt anyway. 

In the present work, no attempt has been made to calculate a contribution to the photon spectrum from 
pre-equilibrium matter. It is, however, unlikely that a long-lasting pre-equilibrium phase is characterized 
by a strong photon emission rate, since strong photon emission indicates frequent interaction processes 
in the medium which in turn would lead to fast equilibration. Furthermore, it is plausible that a pre- 
equilibrium contribution mainly affects the (uncertain) high momentum region of the spectrum, as it is 
characterized by hard momentum scales. Nevertheless, all results discussed in this paper are subject to 
some uncertainties resulting from the poor knowledge of pre-equilibrium dynamics. 

5 Comparison to other works 

Several other scenarios have been investigated by different authors in order to explain the photon spec- 
trum measured by WA98. In |24j . a hydrodynamical evolution model has been used. The favoured 
scenario found in this work uses a very short equilibration time of 0.2 fm/c, corresponding to an (av- 
erage) initial temperature of 335 MeV. The contribution of thermal photons is about 50% of the total 
yield, the rest is prompt photon contribution. 

While the result for the equilibration time is clearly very different from the findings of the present work, 
several other results are similar, among them the weak sensitivity to the phase transition temperature 
Tc (indicating the dominance of early emission phases), the favoured large initial temperature and the 
dominance of the QGP over the hadronic signal. It remains to explain the remarkable difference in the 
conclusions about the equilibration time. 

Note that both in 24 and the present work large initial temperatures Tj > 300 MeV are needed to give 
a sizeable thermal photon emission in the momentum range in question. The choice of To = 0.2 fm/c 
seems mainly driven by the need to create such large Tj. There are, however, two important differences 
between our model and the one in |24| which lead naturally to large initial temperatures for equilibration 
times above 0.5 fm/c in our approach. 

First, we employ an EoS as based on lattice results, which leads to a temperature increase of about 30% 
for a given entropy density as compared to a bag model EoS. Second, the temperatures quoted in 
are obtained using the Bjorken estimate 25 . Our fireball evolution, however, incorporates significant 
longitudinal acceleration of matter, which in essence leads to a peaked initial distribution of energy 
density at central rapidities and hence to significantly larger initial temperatures. 



10 



In |26| . a number of scenarios with different EoS and initial state have been investigated within a 
hydrodynamical description. The reference scenario described there uses a bag model EoS for the QGP 
phase with a transition temperature of 180 MeV and an initial state which leads to a peak initial 
temperature of T™ ax — 325 MeV and an average initial temperature of T(z = 0) = 255 MeV. There is 
clearly a discrepancy between the average initial temperature in |26| and the present work. Due to the 
different space-time expansion patterns of the hot matter in the (averaged) present calculation and a 
hydrodynamics evolution, this issue could possibly best be clarified by comparing the amount of four- 
volume corresponding to a given temperature instead of comparing the average at a given r. At the 
moment however, this has to be regarded as an open question. 

In addition, there is a larger contribution to the photon spectrum of matter with temperatures below 
200 MeV observed in (about 40% of the total yield) as compared to our model. This is presumably 
caused the different choice of the EoS, which in the bag model case leads to faster cooling and to a 
long-lasting mixed phase, in essence reducing the weight of contributions from high temperatures and 
enhancing the low- temperature yield. 



6 Summary 

A measurement of hard real photon emission provides a good opportunity to study the early evolution of 
a fireball. Due to the temperature dependence of the photon emission rate, the QGP phase is expected to 
dominate over the contribution from the hadronic phase. We have used the leading order photon emission 
rate from the QGP, along with estimates of the impact of our phenomenological quasiparticle picture 
on this rate, in a simple model for the fireball evolution to calculate the resulting photon spectrum. As 
this model was fixed beforehand, we have not introduced any new free parameters. 

For the 'standard choice' of the equilibration time To, the model has been able to give a good description 
of the data. Nevertheless, we have tried variations of this quantity in order to work out constraints. 
With the help of an estimate for the contribution of photons from initial, hard processes, we found for 
the equilibration time 0.5 fm/c < tq < 3 fm/c. This would imply that a QGP phase must be present in 
the evolution, at least within the present model. 

Overall, our picture of the spacetime evolution of the fireball finds now support from both the low 
momentum (late time) and the high momentum (early time) region of the evolution. More precise 
future measurements and calculations can be expected to tighten the constraints on the equilibration 
time. 
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